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Abstract 

The general principles of the choice of the reduced description parameters of 
nonequilibrium states J a (t) an d the construction of the nonequilibrium statistical 
operator (NSO) p(t) are discussed. On the basis of Kavasaki - Gunton projection 
operator an explicit non-Markov form of NSO with respect to j a (t) parameters is 
obtained and the equivalence of the method to the other methods using boundary 
conditions for p(t) or "mixing" operator is shown. An exact NSO transformation to 
Markov form with respect to 7«(i) at arbitrary dependence of Vt perturbation on 
time is made. The generalized kinetic equations not containing the memory with 
respect to j a (t) are obtained being convenient for practical application. 
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1. Introduction 



The existing method of theoretical study of irreversible processes [1-3] are based on 
N.N.Bogoliubov concept [1] of sequential reductions of parameter number necessary for 
description of nonequilibrium states of macroscopic system in the process of its evolution 
in time. In accordance with this concept a quasiequilibrium states is determined in the 
system at t ^> r (t - is the chaotization time) that can be characterized by giving 
average values j a (t) = Sp{p(£)7a} varying comparatively slowly in time and generally 
speaking in space. {7»} is a set of limited linearly independent operators (a is the complex 
of discrete and continuous indices on which 7 a operators depend). 7 Q (t) is generally 
called the generalized thermodynamical coordinates (GTC) of system. Nonequilibrium 
statistical operator (NSO) p(t) giving a rough description of the system at t 3> t is, in 
the general case, a functional depending on operators 7 a and GTC 7 Q (t). 7 a (t) functions 
are determined by solution of the closed system of generalized kinetic equations (GKE) 
obtained by averaging of operator equations of motion for 7 Q over NSO p(t). 

The choice of 7 a operators should be made individually for each specific system taking 
into account the hierarchy of the interactions existing in it. 

2. Choice of GTC 

Let us consider space-uniform nonequilibrium system with time-dependent Hamilto- 
nian H t : 

H t = H + V t ; V t = V + H F (t) < H (2.1) 

H — is the main Hamiltonian defining a set of operators {7q} and V t contains relatively 
weak interaction including those with variable external fields (addend Hpit)) that do not 
affect the choice of operators 7 a . 

Generally, the degree of freedom of macroscopic system can be divided into separate 
quasi-independent groups called subsystems. Then, H is the Hamiltonian of free sub- 
systems including as a independent subsystems the effective multipart interaction H e ff 
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responsible for chaotization time To; V is the interaction between subsystems. Such situ- 
ation is realized in the problems of magnetic resonance in solids [4]. 

At investigation of the kinetics of condensed medium on the basis of the idea of 
quasiparticles [5], as well as in gases [1] H is the Hamiltonian of free quasiparticles and 
the interaction between them is described by operator V. In this case tq is of the order 
of collision time and does not depend neither on H , nor on V value. 

The main requirement for 7a, operators is [3]: 

iLola = i [la, H ] = % OiapifS (2.2) 

in 

a a p is the matrix of c numbers defined by properties of H operator symmetry. For solid 
spin systems the following requirement should be added to (2.2) 

[7a, H eff ]=0. (2.3) 

The idea of (2.2) and (2.3) conditions is that the set of {7 a } is formed by the main 
Hamiltonian H independent of V t perturbation. Operators 7 a while commutating with H 
are to be closed on themselves, i.e. they should commutate with the effective interaction 
when it is included into Hq, but not necessarily commutate between themselves. 

According to (2.1) - (2.3) the evolution of macroscopic parameters j a (t) in time in- 
cludes the dynamic change at the expense of H (being considered exactly) and H F (t) <C 
H (being considered according to the perturbation theory), as well as the relaxation 
change at the expense of V <C H with relaxation time r r , and in accordance with the 
principles of reduced description r r ^> t . 

3. The method of construction of nonequilibrium statistical operator 

The simplest method of NSO construction is the application of Liuville equation with 
infinitively small source [6]: 

+ iL(t)}p(t) = -e{p{t) - P(t)p(t)}; e - +0; (3.1) 
where Liuville operators are determined by the expression: 
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L(t) = Lo + L v (t); 
iL A = -^\A,Ho]; 
iL v (t)A = ^[A,V t ]. 

Projection operator P(t) determines the degree of rough description of nonequilibrium 
states of the system and the structure of NSO pit) at t ^> r . It is convenient to use 
Kavasaki - Gunton projection operator as a P(t) [7]: 

P(t)A = Pq {t)SpA + £ p^{S P (i a A) - 7a(t)SpA}; (3.2) 

where 

a 

Q q = Spexp{-J2 F a(t)ia}; (3.3) 

a 

p q (i) is the quasiequilibrium statistical operator depending on t implicitly through GTC 
7 a (£). The functional connection F a (t) = F a {j a (t)) is determined by matching conditions 

la(t) = Sp{p{t)i a } = Sp{p q (t)i a }, (3.4) 

providing the coincidence of true and quasiequilibrium values of macroscopic parameters 
and the conjugation of two sets of j a (t) and F a (t) functions in the sense of nonequilibrium 
thermodynamics [2]. 

According to (3.2) and (3.4) P(t) has the following properties: 

P(t)p{t) = p q (t); 

[ ) dt - dt ' 

P{t')P{t)A = P(t')A. (3.5) 
Eq. (3.1) is similar to usual Liuville equation: 

{^ + iL(t)}p(t) = (3.6) 



and to the boundary condition of delay type [6,8] 

/t+T ft+T 
L(t')dt'}p(t + r) — w_oo Texp{i J L(t')dt'}p q (t + r). (3.7) 

Thus, the introduction of infinitively small source into Liuville equation is the conve- 
nient method of selection of delayed solution of this equation that depends on 7 a (t) and 
satisfies the boundary condition in (3.7). 

In the Appendix the equivalence of the boundary condition (3.7) and the mixing 
operations are shown: 

exp(iL T)p(t) — > T ^-co exp(iL T)p g (^(t)), (3.8) 

that is used in the method [3,9]. 

This fact is very remarkable the more so that we have the significant difference in 
conditions (3.7) and (3.8). If in (3.7) evolution of the system is considered from moment 
t to moment t + r along the phase trajectory, i.e. with the whole Hamiltonian H t during 
which are changed both the dynamical variables (multiplier Texp{i // +r L(t'dt'}), and 
the macroscopic parameters (7«(t) — > 7 Q (t + r)), in (3.7) the mixing operation describing 
a free evolution of the system takes place at fixed values of 7o,(t). 

The solution of Eq. (3.1) has the following form [6,8]: 



p{t) = p q (t) - J° dre eT Texp{i L(t')dt'} ( ^Llll ^ + lL {t + T )p q {t + r). (3.9) 

The use of explicit form of NSO (3.9) in the Appendix is difficult because of its complex 
structure. In this respect the integral equation for p{t) [10], derived in the absence of 
variable fields for 7 Q operators satisfying (2.2) is much easier: 

pit) = Pq (t) - f dre^ L ^{l - P R (t + r))iL v (t + r)p(t + r), (3.10) 



where 
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is the Robertson projection operator [11]. According to (3.2), P R (t)A = P(t)A, if 
SpA = 0. 

At e = Eq. (3.10) coincides with the equation obtained by the other method [12] 
when studying the system evolution under the influence of external fields variables assum- 
ing that at r — > — oo the equilibrium condition of the system was disturbed by adiabatic 
application of variable field. 

In (3.9) and (3.10) p(t) is non-Markov functional of 7 a (t) as it depends on all values of 
la{t') t' < t. As a consequence, the non-Markov integral-differential generalized kinetic 
equations are obtained even in the absence of high-frequency fields (H F (t) = 0) [2]: 

= SpiJ^- = S PP {t)iL{t)i a = L a (t), (3.11) 

where L a (t) is the collision integral. 

On the other hand, in accordance with the requirement of the completeness of the 
reduced description of irreversible processes a nonequilibrium condition and hence pit) at 
some moment of t (from macroscopic point of view) i.e. for t ^> r , should be determined 
by giving 7 a (t) parameters at the same moment of time. In particular problems the 
memory seems to be insignificant due to the use of operator with (2.2) and (2.3) properties 
as 7q, operators and to the weak interactions V t [4]. In stationary case V t = V ^ f(t), 
as it is shown in [13], one can, in principle, obtain Markov presentation of NSO in any 
finite order with respect to V by reconstruction of a number of perturbation theories with 
respect to weak interaction V [10]. 

Taking into account that in widely and successfully used methods of GKE construc- 
tion the Markov character of NSO pit) dependence on 7 a (t) is supposed from the very 
beginning [1,3] it is more natural to formulate the Markov form of NSO II method [6] in 
general form without using a definite order of perturbation theory for V t . The limitation 
of the class of Eq. (3.1) solutions to Markov type eliminates the difficulties in solution of 
non-Markov equations (3.10) in high orders of perturbation theory with respect to V t . 
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4. Conversion of NSO to Markov form 



Let us introduce the operator of shift in time 

Q oo Qn 

exp{r—Mt) = g ^^(*) = ^ + T) - 
Then Eq. (3.10) has the following form: 

Pit) = Pq (t) - f dre^ L ^e^{(l - P R (t))iL v (t)p(t)}. (4.1) 

J — oo 

Let us proceed from the superposition 

p(t) = p(j{t),t); t>r , (4.2) 

assuming that p(t) at t ^> r depends on t not only through 7 Q (t) value at the same 
moment of t, but by V t . 

Note, that the time dependence of NSO just of such kind is obtained as a result of 
direct reconstruction of a number of perturbation theories [13], thus the supposition (4.2) 
can be considered as a mathematical formulation of the possibility to construct Markov 
with respect to j a (t) NSO expansion into series with respect to V t perturbation. In this 
case the effect of d/dt operator on arbitrary Markov type f(^(t),t) functional of 7 Q (£) is 
determined by the relation: 



|/W f ). f )=(E^^ + |)/W').^WTW.<) + |)/(TW,<), («) 

where means the partial derivative with respect to real time included in f(^(t),t) 
functional. 

Differentiating (3.4) with respect to time and using (3.1), (2.2) and (3.5), we obtain 
GKE 



d la {t) 



i 



dt p 



^a Q/3 7^(t) + 
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Sp(p(j(t),t)iL v i a ) = L a (j(t),t), (4.4) 

where L a is the collision integral of Markov type with respect to 7a (i). 
According to (4.3) we can write 



D(l(t),t) = E L «(7^)/- = A>(t) + A,(7,*)l 

where 

d 

D v = Y,Sp{p{l,t)iL vla )^-. (4.5) 

a Ola 

Introducing the symbols: 



b = pin,t) -p,(t); 

A=-(l-P R (>y))iL v (t)p(>y,t); 
H= -(e + iL + D (j)); 

U=-(D v (^t) + ^- (4.6) 
% at 

in view of Eqs. (4.2), (4.3), (4.5) and the rearrangement ability of (e + iL ) and D("f,t), 
operators, Eq. (4.1) can be written as: 



B = [° dre* H+u ^A. (4.7) 
Eq. (4.7) is equivalent to integral equation [3] 



B 



= [ dre lHr {A-iUB). (4.8) 

J — oo 



In previous notations, taking into account that exp{rD {^)) operator describes free 
evolution of (V t = 0) from t moment to t + r moment at which according to (4.4) 7«(t) — > 
7a (t + r) = (e iar 7(t)) a = e rA,(7) 7 a (t) Eq. (4.8) takes the following form: 



fat) = p q (t) - f dre^ L ^{iL v {t)p{ 1 ,t) + 

J —oo 
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+ E ^^p(p( 7 , t)iL v {t) la ) + ^p( 7 , t)} 7 ^ 7 , (4.9) 
where the equality 

A,(7,*)p,(7) = E%^^{P(7)^W7J = -P R ( 7 )^Wp(7,t). 

is used. 

Similarly, instead of (4.6) we use the notations (symbols): 



b = pin,t) - p q in)\ 

A=-(l-P R (>y))iL v (t)p(>y,ty, 

H = -{t + iL, + D {) { 1 ) + ^-)- 
% at 

U=-D v ( 1 ,t), (4.10) 
% 

we obtain: 



p(7,t)=p 9 ( 7 )- f dre^ +iL ^{iL v (t + T)p( r t + r)+ 

J — oo 

+ E 9p(7 , t + r) ^p( 7 , t + r)il v (t + r)f Q } 7 ^ 7 . (4.11) 

It is easily seen that Exps. (4.9) and (4.11) are strictly equivalent and automatically 
satisfy the matching conditions (3.4). 

When e = 0, (4.9) coincides with the equation obtained directly from Liuville equations 
using boundary conditions (3.8) in low-frequency region ujtq <C 1 of variable field [3], and 
Eq. (4.11) coincides with the equation derived by similar method for arbitrary values of 
lut parameter [14]. In stationary condition (H F (t) = 0), when NSO p(t) does not depend 
explicitly on time at e = Eqs. (4.9) and (4.11) coincides with the basic equation of the 
method [9]. The same result assuming p(t) = p( 7 (£)) at t 3> t can be obtained directly 
from Eq. (3.1) [15]. 

In our method Eqs. (3.10), (4.9) and (4.11) are similarly applicable at arbitrary 
values of u t parameter. Nevertheless, in the low-frequency region when ut <C 1 and 
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correspondingly p \£ ' ~ ut term is small it is convenient to use Eq. (4.9) and at high 
frequencies ut > 1 is preferable to use Eq. (4.11). 

In general case at arbitrary lut GKE (4.4) are differential with respect to t, i.e. belong 
to Markov type, though they contain the memory with of respect to variable field. At 
ujt <C 1 the memory with respect to field disappears and the evolution of the system in 
time in the variable field is strictly Markov type. 

5. Discussion of results 

Let us discuss the conditions of applicability of the obtained equations. According to 
(3.6) p{— oo) = p q (—oo) = p (p - is the equilibrium statistical operator of the system). 
To form the quasiequilibrium state ~ r time is required. Therefore NSO p(j, t) and Eqs. 
(3.10), (4.9), (4.11) describe the system at t > r . 

Let us assume that in remote part the system was in equilibrium (p(— oo) = p ) an d 
then the equilibrium was disturbed by adiabatic application of variable fields. As is shown 
in [3,12] at such boundary condition the limitation of t 3> t is removed and Eq. (3.10) 
and hence Eqs. (4.9) and (4.11) following from it describe the system adequately all t 
time. 

As for the frequencies u and amplitudes lui of variable fields Eqs. (3. 10), (4.9) and (4.11) 
are applicable at arbitrary frequencies and at the amplitudes for which H F (t) <C H Q . 

Macroscopic parameters 7 Q (t) usually are introduced for description of nonequilibrium 
states of systems, the Hamiltonian of which does not depend on time. At the action of 
variable fields the time-dependence of Hamiltonian is eliminated by the transition into 
rotating coordinate system [16] or by considering the fields from quantum point of view 
[17] and only after this 7 Q (t) is introduced. When 7 a operators commutate with H i.e. 
a a/ 3 = such approach preserves the common thermodynamic sense of 7«(t) values and 
is the most consistent. But such approach is not always applicable, therefore one should 
discuss the problem of possibility to use ^ a (t) parameters in general case. 

First, let us assume that a a p = 0. In the low- frequency region r y a (t) parameters, being 
the slowly changing time functions preserve the common thermodynamic sense. At each 
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instant of t time the system is "adjusted" to the instant value of Hamiltonian H t (2.1). 
Then the whole difference from the stationary case is reduced only to that for t ^> tq 
the statistical operator p(t) depends on t not only through 7 Q (i), but explicitly through 
the variable field F(t) and its derivatives F(t), F(t), . . . [3]. As it mentioned in [3], for 
ujtq <C 1 the kinetic processes are in fact the Markov type and (4.9) is obtained from (3.10) 
by direct expansion of V t+T and 7 a (t + r) values into power series r. In our aforementioned 
approach, Eqs.(3.10) and (4.9) are strictly equivalent for any values of ut parameter. 

In the high-frequency region ujtq > 1 the non-Markov character of kinetic processes, 
taking place in the system under the action of variable fields will manifest itself. In this 
case the expansion in terms of ujtq parameter seems impossible. 7«(t) parameters have 
no longer the direct thermodynamic sense as they contain the rapid, but low-amplitude 
(because H F {t) <C H ) dynamic change induced by high-frequency field. Nevertheless, 
7 Q (t) parameters can be still used for description of quasiequilibrium states, as, according 
to (3.2), owing to the "mixing" action of H , the system will have time to be "adjusted" 
to the instant nonequilibrium values of 7 Q (t) parameters. 

When a a p ^ and V t = 0, 7 Q (t) values in accordance with (4.4) are changed by 
the law 7(t) = e w *7(0) with arbitrary frequencies determined by eigenvalues of (e tat ) afj 
matrics. Therefore, j a (t) are not "slow" thermodynamic values even in the absence of 
variable fields. However, it is remarkable that in the general case a a p ^ 0,ut > 1 the 
"adjusting" of the system to the instant values 7 Q (t) takes place and the dependence of 
p(^(t),t) on 7a(0) as we U as Eq. (4.11), are of Markov character with respect to 7 a (t), 
though they contain the memory with respect to variable field. 

Conclusion 

The carried out study shows the full equivalence of the boundary condition (3.7) and 
mixing operation (3.8). Thus, we prove the strict equivalence of the method [3,9,12] to 
the variety of NSO method [6] based on Liuville equation with infinitively small source 
(3.1) and 7 Q operators with (2.2) properties. 

The main conclusion of the work is that the method NSO II [6] can be formulated in 
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Markov form with respect to j a (t) parameters even in the presence of high-frequency fields. 
As a result, the use of high order perturbation theory with respect to small interactions 
in the system is considerably facilitated. 
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Appendix 

Let us show the equivalence of boundary condition (3.7) and mixing operation (3.8). 
Liuville equation (3.5) can be written as 



|- + iL )5p(t) = -(1 - P R (t))iL v (t)p(t); 5p = p- Pq . (Al) 



Let us multiply (Al) by integrating factor exp(iL t) and integrate it from — oo to t. 
Note, that the solution of Eq. (3.1) has the [6,8]: 

pit) = p q (t) - J° dre"T exp L{t')dt'} ^A ? (* + T ) + tL{r)pq{t + r) ^j = 

= p q {t) - f dre eT T(t,t + r)(l - P R (t + r))iL(t + r)p q (t + r), (A2) 

where 

T(t, h) = Texpj- j\l - P R (t'))tL(t f )dt'Y 

and P R (t) is determined in the main text (see (3.10)). 
Using (2.2) and the relation [3,8], 

e- lLoT Pq ( 1 (t))= Pq (e2^ 1 (t)), (A3) 

we obtain: 

e iL ^8p{t) = - f dt'e lLot '(l - P R (t'))iL v (t')p(t'). (AA) 
Taking into account that according to (3.8), 



hm e ^(p(r)-p q (T)) = 0, 
after simple transformation of (A4) we have 

p(t) = p q (t) - f° dre iL » T (l - P R {t + r)iL v (t + r)p(t + r). (A5) 

J —CO 

(A5) coincides with (3.10) at e = 0. The forming factor e eT does not affect the integral 
convergence [13,15] included in Eq. (3.10) and Eqs. (3.6), (3.8) and (3.6), (3.7) are fully 
equivalent as they result in coinciding integral equations (3.10) and (A5). 
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